Many-flavor Phase Diagram of the (2 + l)d Gross-Neveu Model at Finite Temperature 
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We study the phase diagram of the Gross-Neveu model in d — 2 + 1 space-time dimensions in 
the plane spanned by temperature and the number of massless fermion flavors. We use a functional 
renormalization group approach based on a nonperturbative derivative expansion that accounts for 
fermionic as well as composite bosonic fluctuations. We map out the phase boundary separating the 
ordered massive low-temperature phase from the disordered high-temperature phase. The phases 
are separated by a second-order phase transition in the 2d fsing universality class. We determine 
the size of the Ginzburg region and show that it scales to zero for large Nt following a powerlaw, in 
agreement with large- Nf and lattice studies. We also study the regimes of local order above as well 
as the classical regime below the critical temperature. 



I. INTRODUCTION 

The Gross-Neveu model [l| in (2 + l)-dimensions ex- 
emplifies many of the intriguing properties of interacting 
relativistic fermion systems in a planar spacetime. From 
a theory perspective, it is an ideal laboratory to study the 
interplay of initially massless fermions with dynamically 
generated collective excitations. From the viewpoint of 
applications, the study of effective quantum field theories 
of low-dimensional Dirac fermions has become a rapidly 
evolving field in the context of condensed-matter systems 
like graphene 0, topological band an d Mott insu- 
lators [y], as well as the normal-state behavior of high- 
temperature superconductors 0, HI ■ 

The Gross-Neveu model in (2 + l)-dimensions exhibits 
a quantum phase transition to a phase of broken (dis- 
crete) chiral symmetry as a function of the fermionic 
interaction strength. Despite being perturbativcly non- 
renormalizable, it provides a paradigmatic example for 
an asymptotically safe theory Q that is ultraviolet (UV) 
complete due to the existence of an interacting non- 
Gaufiian renormalization group (RG) fixed point 0. 
This fact has been proven to all orders in a 1/Af ex- 
pansion and is supported by lattice simulations 
[T3 - [l6j as well as functional RG studies [l7|- The quan- 
tum phase transition associated with the RG fixed point 
has been studied with a variety of techniques flftl21T ] . 

Asymptotic safety makes the Gross-Neveu model in 
(2 + l)-dimensions a complete theory which is valid on 
all scales. It also ensures a "perfect" predictive power: 
once, a single (scale) parameter is fixed, each and every 
physical observable can be predicted absolutely. This 
statement extends to all dimensions 2 < d < 4, with the 
non-Gaufiian fix ed p oint merging with the Gaufiian fixed 
point for d -> 2 (l7j }. 

From the viewpoint of critical phenomena, the proper- 
ties of the UV fixed-point control the universal features 
of the quantum phase transition. Quantitatively, this 
statement extends to the long-range scaling behavior of 
the chiral order parameter and/or the correlation length 
near the phase transition which arc both governed by the 



critical exponents of the RG fixed point. Also for these 
critical exponents, the quantitative agreement among: the 
various methods is very satisfactory, see e.g. [16l - fl8l |20| . 

Whereas the quantum phase transition in the Gross- 
Neveu model in 2 < d < 4 dimensions is governed by 
the coupling strength as a control parameter, many low- 
dimensional fermion systems can also develop interesting 
features or even phase transitions as a function of the 
flavor number Af, i.e., the number of interacting fermion 
species. For instance, convincing evidence has been col- 
lected for the fact that the (2 + l)-dimensional Thirring 
model features chiral symmetry breaking (xSB) only for 
flavor numbers below a critical value Nf < Nf, CI [28l — 
[3l| . A similar behavior is found for fermion systems 
with gauge interactions such as QED( 2 +i) [32| or even 
QCD( 3+1 ) [33j . The latter class of many-flavor theories is 
currently actively discussed in the context of New Physics 
scenarios beyond the standard-model of particle physics 

In all these cases, non-analytic features as a function of 
the flavor number Nf occur, limiting the validity region 
of the large- Af expansion. The occurrence of a critical 
flavor number in these models arises from a competition 
between symmetry stabilizing and destabilizing fluctu- 
ations. In this respect, the (2 + l)-dimcnsional Gross- 
Neveu model behaves less dramatically as a function of 
Af. Qualitatively, this can be understood from the fact 
that only one relevant bosonic degree of freedom cor- 
responding to a "radial" Z 2 mode cannot counterbal- 
ance the destabilizing fermion fluctuations at strong cou- 
pling. Still, the Gross-Neveu model is not completely 
structureless as a function of flavor number, as is indi- 
cated by several observations: first, in the limit Af — > 0, 
fermion interactions are switched off leaving us with a Z 2 - 
symmetric universality class of (collective) scalar degrees 
of freedom and a corresponding Ising-typc phase tran- 
sition. This implies that the Gross-Neveu universality 
class has to merge with the 3d Ising universality class in 
the limit of small flavor number, approaching a Wilson- 
Fisher fixed point. From the large- Af expansion, it is al- 
ready visible that this merger requires a non-monotonic 
behavior in the universal quantities. Second, the large- 
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N{ limit (by definition) merges exactly with mean-field 
theory which therefore must also hold for the scaling 
of observables at the finite-temperature transition. On 
the other hand, simple dimensional reduction suggests 
that the finite-temperature phase transition is governed 
by order-parameter fluctuations in the 2d Ising universal- 
ity class which is clearly different from mean-field theory. 
This mean-field puzzle has first been studied in [35j us- 
ing lattice simulations, followed by generalized large- Nf 
expansion techniques in [36] . It turns out that its resolu- 
tion requires the existence of singularities in the critical 
exponents in the large- Nf limit. 

In the present work, we study the phase diagram of the 
(2 + l)-dimensional Gross- Neveu model as a function of 
the flavor number and temperature. Because of asymp- 
totic safety, this phase diagram is a pure prediction of the 
theory once an overall scale parameter is fixed. Our em- 
phasis will be on universal aspects of criticality, as well 
as near-critical aspects such as local-order phenomena or 
the width or the critical region. Of particular interest 
are the approach to the large- Nf limit as well as the ex- 
pected non-monotonic features at small flavor numbers. 
All these different aspects can already be studied within 
a comparatively simple nonperturbative approximation 
of the functional RG flow, whence we concentrate in the 
present work on conceptual clarity rather than quantita- 
tive precision. As our approximation scheme corresponds 
to a next-lo-leading order result of a derivative expansion 
of the effective action, higher quantitative precision can 
be expected from systematic inclusions of higher orders. 

Since the physical spin of electrons in condensed- 
matter applications of the Gross-Neveu model typically 
plays the role of flavor, the Nf = 2 model is phenomeno- 
logically most important. For instance, the excitonic 
[37j| or antiferromagnetic instabilities |HH in graphite and 
graphene, respectively, can be associated with quantum 
phase transitions falling into this universality class. Also 
the universal aspects of secondary d- to d + is- wave pair- 
ing transitions in nodal <i-wave superconductors are ex- 
pected to be determined by the Nf = 2 Gross-Neveu 
model. Therefore, special iVf-dependent structures of the 
many-flavor properties which lie near Nf = 2 can also ex- 
ert an influence on physical systems. 

This work is organized as follows. In Sect. HH we in- 
troduce the model and discuss the relevant symmetries. 
Section Mil briefly summarizes the results of a standard 
mean-field analysis, following earlier works in the litera- 
ture, see, e.g., Ref. [12| . Section flV] is devoted to a func- 
tional RG study of the model at next-to-leading order in 
a derivative expansion. The many-flavor phase diagram 
at finite temperature and further results are presented in 
Sect. [V] We conclude in Sect. IVII and list a number of 
technically relevant details needed for the RG analysis in 
the appendix. 



II. GROSS-NEVEU MODEL IN d = 2 + 1 

The classical action of the Gross-Neveu model in d 
space-time dimensions reads 

S$,il>] = / d d x + E Mi-JjfW, \ 

EE J d^j^V+^W) 2 } , (1) 

where Nf determines the number of fermion flavors. The 
fermions interact via a four-fermion interaction term with 
the bare coupling given by g of mass dimension [g] = — 1 . 
In our study of the Gross-Neveu model in d = 2 + 1, 
we shall use a reducible four-component (d 7 = 4) rep- 
resentation for the 7-matrices, as such Dirac spinors ap- 
pear naturally in the condensed-matter applications of 
this model, see Refs. [lj], for details of our conven- 
tions. 

The symmetry which is most relevant for the present 
work is a discrete symmetry with generators given 
by Z2 = {14,75}- The nontrivial transformation of the 
Dirac spinors reads 

i> -> 75^, ip -> -"075, (2) 

implying ipip — > —ipip and thus leaving the action (fTj) 
invariant. In agreement with the literature, we will refer 
to this symmetry as a discrete "crural" symmetry (even 
though it should rather be viewed as a subgroup of the 
axial symmetry familiar from Ad models.) Note that this 
transformation acts simultaneously on all flavors. 

In addition to this discrete chiral symmetry, the (2+1)- 
dimcnsional Gross-Neveu model has a much larger con- 
tinuous flavor symmetry: using the product of Dirac 
matrices 735 := 17375, we can define the projectors 
P35,± '■= (1/2)(1 ± 735)- For each of the corresponding 
spinor components ij)% ± := P35 1 ±ipi, we can define a sepa- 
rate U(iVf)± flavor rotation, such that the action |1| is in- 
variant under a U(-/Vf) + x U(Af)_ continuous symmetry 
group. However, the latter will not play any relevant role 
in the following and remains unbroken in the discussed 
pattern of dynamical symmetry breaking. FRG stud- 
ies of 3d chiral phase transitions with more com plic ated 
symmetry-breaking patterns can be found in [IlMEl- 

Perturbative arguments indicate that the Gross-Neveu 
field theory in d = 4 dimensions is non-renormalizable. 
This implies that the theory has to be defined with a 
physical cutoff, introducing - strictly speaking - infinitely 
many physical parameters that specify the details of the 
physical regularization. In lower dimensions 2 < d < 4, 
the perturbative reasoning is invalidated by the occur- 
rence of a non-Gaufiian fixed point in the RG flow that 
allows to define UV complete renormalized trajectories 
that emanate from this interacting UV fixed point. This 
asymptotic safety is conceptually reminiscent of UV com- 
pleteness due to asymptotic freedom with the only dif- 
ference that the (dimcnsionless) coupling approaches a 
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finite fixed-point value in the UV. A detailed analysis 
\vi\ shows, that a renormalized trajectory is fully fixed 
by fixing an overall scale parameter. In 2 < d < 4, 
this can, for instance, be done by fixing the dimension- 
ed bare coupling g\ at an initial scale A. For suffi- 
ciently strong initial coupling, the system is in the mas- 
sive chiral-symmetry broken (xSB) phase such that the 
scale fixing can also be done by choosing a value for the 
induced fcrmion mass mf. 

In the present work, we will always assume that the 
initial coupling has been chosen super-critically, i. e. the 
system is in the massive %SB phase and all observables 
can be given in units of mf. To compare different theo- 
ries, we keep m.f fixed to the same value for all N{. Inci- 
dentally, the critical coupling is directly connected with 
the non-GauBian fixed-point coupling, implying a stan- 
dard relation between the RG fixed-point structure and 
critical aspects of the chiral quantum phase transition. 

In passing, we note that the non-Gaufiian fixed point of 
the four-fermion coupling approaches the non-interacting 
Gaufiian fixed point in the limit d — > 2. In other words, 
asymptotic safety becomes asymptotic freedom in d = 2. 
Consequently, the phase transition disappears as the sys- 
tem undergoes xSB for any initial finite value of the cou- 
pling with mass generation being associated with "di- 
mensional transmutation" . Even though large parts of 
the analysis can be done in continuous dimensions, we 
concentrate on d = 3 in the following. 



III. MEAN-FIELD ANALYSIS 

Let us start by recapitulating a mean-field analysis of 
the model and its finite-temperature phase transition. 
Such an analysis has already been performed frequently 
in the literature, see, e.g., Ref. fl2j |. We summarize here 
only the essentials as a preparation for the functional RG 
anaysis. If only fermionic fluctuations are considered as 
in the following, the mean-field analysis is identical to 
the large- Nf limit. 

We begin by introducing an auxiliary bosonic 
Hubbard-Stratonovich field a ~ -0-0, leading us to a so- 
called partially bosonized action: 



N { 



m 2 a 2 + ip (i$ + iha) i\> 



(3) 



Since the a field occurs only quadratically, it can be in- 
tegrated out exactly. This leads us back to the fermionic 
action ((T|) if we identify 
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h 2 



in 



2 ' 



(4) 



condensate (i/iip) , it serves as an order parameter for chi- 
ral symmetry breaking and renders the fermions massive. 
From a study of the order parameter for a given Nf we 
can obtain the critical temperature T cr and analyze the 
universal critical behavior close to the phase transition, 
see below. 

As the fermions occur only bilinearly in the partially 
bosonized action ([3]), they can be integrated out, yielding 
a nonlocal bosonic effective action for the Gross-Neveu 
model, 



S[a] = Nf 



- 2 2 

m a 



AT f Trln(^ + iha). 



(5) 



The mean-ficld/largc-Af approximation is defined by ig- 
noring fluctuations of the a field, i.e., Eq. ([5]) defines the 
full theory at mean-field level. As we expect the ground 
state to be homogeneous for all temperatures T and fla- 
vor numbers Nf, we restrict ourselves to a = <jq = const, 
here and in the following. 1 

For constant <xo, the effective action boils down to 
the effective potential, S[a — » cr ] = J U(ao), with the 
ground state satisfying dU (a) / da\ ao = 0. The solution 
is found from the variation of the action ([5]) , yielding the 
gap equation, 



00 



irr 



(2tt) 3 p 2 



h 2 c 



(0) 



Apart from the trivial solution Co = 0, we find a non- 
trivial solution being the global minimum of the effective 
potential for sufficiently large initial coupling g — h 2 /m 2 . 
The precise evaluation of Eq. (|6|) requires, of course, a 
regularization of the linear UV divergence and a corre- 
sponding rcnormalization of the parameters. Considering 
a sharp UV cutoff A 3> hao for simplicity, we adjust the 
initial bare coupling g = h 2 /m 2 at the scale A such that 
the induced fermion mass fhf = hoo is kept constant. 
Ignoring sublcading orders in the cutoff, the renormal- 
ization condition is given by 



mf 



where we can read off the mean-field value of the (non- 
universal) critical coupling for chiral symmetry breaking, 
5cr = 7r2 /(2A) (for a sharp cutoff). The effective poten- 
tial can be obtained by integrating the gap equation. For 
fhf =/= 0, the result is, 





fh 2 \ 




( 7T 2 1\ 






= -A| 








7T 


\ 2Ag) 



U(a) 



do' 7r V 3 2 



1 



(8) 



h 

such that only the ratio 4-^ has a physical meaning at 
mean-field level. The a field in Eq. ^ transforms as 
a i — y — (7 under the discrete chiral transformation. Since 
its expectation value is directly related to the fermion 



1 Inhomogcneous ground states have been found in the GN model 
in d = 2 dimension |42|1 , but only in the presence of a finite 
chemical potential. Indications for similar phenomena in higher 
dimensions have also been collected, see, e.g., |43|445(I . 



4 



where all the scheme dependence is cancelled by the cor- 
responding scheme dependence of the relation between 
the fermion mass fhf and the initial couplings J4(| ■ Fixing 
the physical observable fhf , the regularization scheme de- 
fines a pair (A, 17(A)), such that a different choice for the 
cutoff A is connected with a corresponding adjustment of 
the initial coupling g(A), such that fhf is unaffected. 

At finite temperature, the same analysis in the Mat- 
subara formalism requires to replace the momenta p — > 
(u! n ,p), with the fermionic Matsubara frequencies u) n = 
(2n + 1)itT, and the spatial momenta p. Performing the 
Matsubara sum, the analogue of the gap equation © (for 
(T ^0) reads BIO, 



1 = 45 



d 2 p 



(27T)2 ^2 



1 - 2n F [ sip 1 + h 2 a 2 ) 



(9) 

where tif{z) = l/(c ,3z + 1) denotes the Fermi distribu- 
tion function at inverse temperature j3 = 1/T. For the 
system being in the chirally broken phase with fermion 
mass fhf > at zero temperature, the gap equation can 
be reexpressed in terms of the finite-temperature fermion 
mass, 



TO£ t = m{ — 2T In ( 1 



(10) 



For fixed zero-temperature fermion mass fhf, this tran- 
scendental equation can easily be solved numerically for 
a given temperature, yielding a thermal fermion mass 
m^T which is monotonically decreasing with temper- 
ature. The mean-field critical temperature T cr where 
fhf X—T CT = can analytically be read off from Eq. (fTU)) , 



rrif 
21n2 



0.7213 fhf. 



(11) 



Analogously to Eq. ([8]) , the thermodynamic effective po- 
tential can be computed by integrating the gap equation 
in terms of poly logarithms, yielding [4y, |47j 



2h\a\T 



+ ^-Li 2 (-e-Rkl/r) 



(12) 



0713 



where the last term involving the Riemann C function has 
been chosen such that U{a = 0; T) = 0. It is instructive 
to expand the effective potential in powers of the mean 
field: 



U(a;T) = 



h 2 In 2 



(T - r cr ) a 2 



16ttT 



576ttT 6 



(13) 



Since only even powers of a appear, the effective potential 
is analytic in the Z 2 invariant a 2 at finite temperature 



- in contrast to the T = solution, where we observe a 
(ct 2 ) 3/2 dependence, cf. Eq. (8). 

Finally, we can use these findings to extract thermo- 
dynamic critical exponents in the mean-field approxima- 
tion. The exponents a and (3 follow straightforwardly 
from the gap equation © and the effective potential 
Eq. (T2), Eq. (13). Expanding Eq. (TO) close to T cr to 
second order in mfj, we obtain the relation [l2T | 

m f , T 2y/m^y/T CI - T, (14) 

from which we can read off the scaling exponent of 
the order parameter f3 = 1/2. Plugging the result for 
the temperature-dependent ground state into the effec- 
tive potential, we obtain the free energy f(T) per flavor 
species and per unit area. Expanding the thermal contri- 
bution to fourth order in fht t T and differentiating twice 
with respect to temperature, the specific heat per unit 
area can be found as 11211 



d 2 



8(ln2)^ 



T 



T=T C , 



(15) 



Since the free energy vanishes in the symmetric phase, 
this implies a = 0, i.e., a jump in the specific heat across 
the transition. For the exponents S and 7, the analy- 
sis has to be repeated with the inclusion of an exter- 
nal source, whereas the correlation exponent v and the 
anomalous dimension rj a can be derived from the large- 
Nf limit of the induced bosonic propagator. For instance, 
the inverse correlation length, corresponding to the in- 
duced bosonic mass, £ -1 = m a , turns out to be propor- 
tional to the fermion mass, the temperature-scaling of 
which has been worked out in Eq. (fT4|) above. The corre- 
lation exponent £ ~ |T cr — T\~ v thus yields the standard 
mean-field value v = 1/2. As it should be, also the other 
exponents yield mean-field values, such that the critical 
behavior of the Gross-Neveu model in (2 + 1) dimensions 
in the large- Nf limit can be summarized by 

a = 0,/3 = ^ 7=1, 5 = 3, v = X -, r? a = 0. (16) 

For further details, we refer to [l2|. These exponents 
do not satisfy hyperscaling relations, as can be expected 
from the fact that bosonic fluctuations which establish 
such relations are completely suppressed in the large- Nf 
limit. However, this still raises a question concerning 
the analyticity of the large- Nf limit: at any finite flavor 
number, universality suggests that the critical behavior 
of the thermal phase transition should be governed by 
order-parameter fluctuations. Since the theory has a Z 2 
symmetry, we expect the phase transition to be in the 
Ising universality class. This is not reflected by the large- 
Nf limit given by Eq. (|16[) . Concerning the critical expo- 
nents for the transition to the high-temperature phase, 
the approach to the large- Nf limit must therefore be non- 
analytic in Nf. 

Actually in all of the above-listed mean-field results, 
the dependence on the flavor number Nf naturally 
dropped out. Non-trivial flavor information thus requires 
the inclusion of both fermionic and bosonic fluctuations. 
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IV. FUNCTIONAL RENORMALIZATION 
GROUP ANALYSIS 

For our nonperturbative study of the many-flavor 
phase diagram of the Gross-Ncveu theory, we employ 
a functional RG equation for the quantum effective ac- 
tion. For this, we consider a scale-dependent effective ac- 
tion Tk , governing the field dynamics at an infrared (IR) 
scale k. The so-called effective average action Tk already 
comprises all quantum fluctuations from the UV cutoff 
A down to k. In the limit k — J- 0, the effective average 
action is identical to the full quantum effective action, 
i.e., the 1PI generating functional, whereas it approaches 
the classical action in the limit k — > A. The k evolution 
of Tfc is determined by the Wetterich equation j48| . 

d t r k m = I S Tr{[ri 2) [$] + R k ]- 1 (d t R k )}, d t = fcA 

(17) 

(2) 

The quantity T y k ' denotes the second functional deriva- 
tive with respect to the field variable $, collecting all 
fcrmionic and bosonic fields. The momentum-dependent 
function R k is a regulator, suppressing IR modes below a 
scale k. For reviews of the functional RG adapted to the 
present context, see Refs. j~H)l : "~)">i . At finite temperature, 
the effective action is directly related to thermodynamic 
quantities. For instance, the corresponding effective po- 
tential at its minimum corresponds to minus the pressure. 
In the present work, we intend to study the many-flavor 
phase diagram in a nonperturbative fashion using our 
continuum RG approach. For this, we use as an ansatz 
for the effective action: 

+N i U(o 2 )}, (18) 
I 



where U denotes the effective potential for the order pa- 
rameter. The couplings and the wave-function renormal- 
izations Z a and are assumed to be scale dependent. 
This ansatz is related to a systematic derivative expan- 
sion at next-to-leading order and has already proved to 
quantitatively describe the quantum-phase transition at 
zero temperature [13, HE, . For simplicity, the present 
ansatz ignores possible anisotropic wave function renor- 
malizations due to the presence of the heat bath, see, 
e.g., 

At the cutoff scale k = A, the initial conditions for the 
couplings are fixed such that the action ^ is recovered 
at this scale. This is ensured by the choice 

Z^A < 1, Z^a = 1, h\ = g\m\, U A (a 2 ) = \^m\a 2 . 

(19) 

From our ansatz (fT5)) . we obtain the following flow 
equation for the dimensionless order-parameter poten- 
tial u{p) = k- d U{a 2 ): 

dtu = —du+(d—2 + 7] a )u'p—2d 1 VdlQ F " >d (2h 2 p;r]^) 
+^-2v d li(u'+2pu";Va), (20) 

where 

p = \z a k 2 - d a 2 , (21) 

and the dimensionless renormalizcd Yukawa coupling 

h 2 = Z^Z-'k^h 2 . (22) 
The corresponding flow for the Yukawa coupling reads, 



d t h 2 = {d - 4 + 2r^ + Va )h 2 + —8v d h 4 l[ F 1 B)d {2h 2 K, [u'(k) + 2ku"(k)} ; m „ Va ) 

JVf 

— - U%ku"{k) + 32kV"(k)) v d h* l[ F 2 B)d (2h 2 n, [u'(k) + 2ku"(k)] ; r^, r/ CT ) 

_ J_32 Wd h 6 k l 2 FB)d (2h 2 Kl [u'(k) + 2ku"{k)] ; W , 77,), (23) 
iVf 



where the right-hand side is understood to be evaluated at the fc-dependent minimum of the effective potential, 
i.e., p m in = k = in the symmetric and p m i n = n =/= in the symmetry-broken regime. Here, k is related to the 
dimensionless renormalizcd expectation value of the scalar field, 



K - 1Z ^ (24) 
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satisfying u'(k) 



in the broken regime. The same conventions apply to the anomalous dimensions 

1 c, V d, 



—8^h 2 m[ F 2 B)d {2h 2 K, [u'{k) + 2ku"(k)} ; n4> , Va ) . 
iVf a 

a 

— 8— 

Nf d 



mf (2h 2 K; r]^) — 2h 2 nm^ (2h 2 k; 77^) 



k (3u"{k) + 2ku'"{k)) 2 mf ([«'(«) + 2ku"{k)] , 0; 7 h ) . 



(25) 



(26) 



All so-called threshold functions I, m occurring in these 
flow equations depend on the details of the regulator. 
Formulas for general regulators as well as explicit ex- 
pressions for the linear regulator [57| used in this work 
are given in the appendix. These threshold functions 
describe the behavior of the regularized 1PI diagrams 
across mass thresholds and thus describe the decoupling 
of massive modes from the flow. They also carry the 
temperature dependence and thus encode also the decou- 
pling of higher Matsubara frequencies as a function of the 
dimensionless temperature r = 2nT/k. A diagrammatic 
representation of the flow equations corresponding to our 
ansatz for the quantum effective action is given in Fig. [T] 



V. MANY-FLAVOR PHASE DIAGRAM AND 
CRITICAL BEHAVIOR 

In order to keep the present analysis simple, we solve 
the flow of the potential in a straightforward polynomial 
expansion, 



71— 1 



and u(p) 



00 * 
—r(P- 



(27) 



in the symmetric regime and broken regime (k > 0), 
respectively. Moreover, we ignore higher order couplings 
with n > 3. In fact, higher-order approximations, though 
widely used in the literature, appear to induce artificial 
instabilities in the present model in a (physically less 
relevant) temperature interval below T CT (Nf). We have 
checked that this approximation does not oversimplify 
the model by verifying all known limits. For instance, as- 
pects of the large- TVf/mean- field limit such as the critical 
temperature are reproduced by our simple approximation 
with an error of ~ 3% (T cr ~ 0.7439 TOf in comparison 
with the exact result T cr ~ 0.7213 m f , cf. Eq. (fTTjl). We 
assume that this systematic error is also indicative for 
the error size at finite Nf. 

The following results are obtained from numerically in- 
tegrating the RG flow in the plane of temperatures and 
flavor numbers (T,Nf). Apart from the errors of the nu- 
merical routines which are much smaller than the above- 
mentioned systematic error, we also have to make sure 
that the initial conditions of the flow lie on the "renormal- 
ized trajectory" , i.e., the line of constant physics (A, g(A), 



dtu oc —du + (d — 2 + rj a )u' p 
+ 3H 



\ / '0 



d t h 2 oc (d - 4 + 2jty + r h )h 2 
1 
1 




'■0 



(F)d 




Va OC d t 



j(FB)d ,(FB)d 
'1,1 i '2,1 



0, 



,(FB)d 
'1,2 



(FB)d 




(F)d {F)d 

mo . m\ 



N r d t - -< D>- 



mi 



Figure 1: Diagrammatic representation of the flow equations 
Eq. ((21]), Eq. (J23J) , Eq. ^ and Eq. (25]). The 1PI dia- 
grams represent quantum fluctuations renormalizing the ver- 
tex functions at the RG scale k. Full (dashed) lines repre- 
sent our ansatz for the fermionic (bosonic) propagator modi- 
fied by a regulator term Rk- Black disks denote the Yukawa 
vertex, while gray triangles denote vertices ~ U"'(a) which 
contribute only in the symmetry-broken regime of the flow. 
Below the 1PI diagrams, the labels (I, m) of the associated 
threshold functions are given, see appendix for details. The 
<9t-derivative acts on the regulator dependence of the regu- 
larized propagators and yields regulator insertions dkRk, as 
denoted by crossed circles in the flow-diagrams for the effec- 
tive potential. Note that all diagrams with an internal bosonic 
line are suppressed by factors of 1/JVf. 



that relates initial conditions for different choices of the 
cutoff A to one and the same physical system in the long- 
range limit. Moreover, this trajectory must also be on 
top of or at least sufficiently close to the critical hyper- 
surface emanating from the non-Gaufiian fixed point in 
order to interconnect the long-range limit with the UV 
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Figure 2: Phase boundary of the Gross-Neveu many-flavor 
phase diagram as within our ansatz as obtained from set I of 
initial conditions. We performed temperature scans for flavors 
N { = 1, 2, 3, . . . , 19 and N f = 20, 30, ... , 100. Also shown are 
fits to the phase boundary in the large- Nf and small- Nf region, 
see Eq. and Eq. poj) . 



complete theory. 2 In practice, we keep the UV cutoff 
fixed at A = 100 a. u. (auxiliary units). Then we choose 
two different sets of initial parameters for the Yukawa 
coupling and the mass parameter that are close to the 
critical manifold. Integrating the flow down towards the 
IR and expressing all dimensionful quantities in terms of 
the same scale, such as the rcnormalizcd fcrmion mass, 



(28) 



has to lead to quantitatively identical results. 3 The first 
set (set I) of parameters consists of choosing Ha = 18.0 
and tuning the bare bosonic mass term A2,a such that 
m { = 1 a.u. In the second set (set II), we choose = 6.0 
and tune the bare boson mass such that mf = 5 a.u. The 
initial bosonic wave function renormalization is set to 
Z a a = 1.0 x 10~ 3 , such that the bosonic field is (almost) 
purely auxiliary at A and the corresponding four-fcrmion 
vertex is (almost) pointlike. Any difference between the 
results for these two different sets are due to nonuniver- 
sal finite cutoff effects, the influence of irrelevant opera- 
tors and the corresponding limitations of the truncation. 
Again, the differences which we find are smaller than the 
above-mentioned systematic errors. 

A typical flow at T = starts at k = A in the symmet- 
ric regime, and is then driven by fcrmion fluctuations into 



Whereas chosing a trajectory exactly on the critical hypersur- 
face is essential for a UV-complete limit, it is numerically not 
relevant for the long-range physics, as trajectories in the attrac- 
tive domain of the fixed point are attracted towards the critical 
hypersurfacc exponentially fast. 

Note that at the mean-field level Zj, = 1, rendering the renor- 
malization factor trivial. In that case the notions of rcnormalizcd 
and unrcnormalized mass coincide. 



the broken regime where the flow freezes out, because of 
mass generation. At finite Nf, the fermion fluctuations 
are somewhat counterbalanced by the boson fluctuations. 

Loosely speaking, finite temperature T > favors 
bosonic fluctuations. This is a generic mechanism as the 
fermions acquire non-zero thermal masses because of the 
anti-periodic boundary conditions, whereas this thermal 
decoupling does not apply to the zeroth Matsubara mode 
of the bosons. For < T < T CI , the zero-temperature 
condensate generated by the fermions becomes somewhat 
reduced by the bosonic zeroth Matsubara mode. In this 
picture, T cr corresponds to the temperature where the 
bosonic zeroth Matsubara mode is capable of driving the 
scalar condensate to zero again such that the flow ends 
again in the symmetric regime. From this point of view, 
it is also clear that the flow even enters the broken regime 
at intermediate scales for T > T CI , but, of course, ends 
up in the symmetric regime. This phenomenon of lo- 
cal order is also generic for Yukawa-like systems both in 
relativistic (5(| as well as in non-relativistic systems [58| . 

Our results for the phase diagram in the (T, Nf) plane 
up to Nf = 100 are shown in Fig. [2] Also shown is a fit 
to the large- Nf behavior (dashed line) 



T cr (JV f ) 



0.72686- 



rrif 



1.55098—, 



Here, we did not fix the leading 
exact mean-field value, see Eq. 
order (NLO) term ~ 1/Nf is 
with large-iVf techniques [36 
is known. The NNLO term has 
as (In N { )/Nf, sec Ref. [Hj. 

For practical purposes, we also present a fit to the 
smaller Nf region (solid line), 



(fit:JV f = 20,..., 100). 

(29) 

ig order term (A^ ) to the 
(flT|) . The next-to- leading 
in principle computable 
though no explicit value 
been determined to scale 



T CI (N ( ) 



mf 



0.695994- 



1.04127 



1.11289 2.18236 2.43492 



Nf 

(fit: Nf = 2 



Nf 



Nf 



19). 



(30) 



In Fig. [3] as well as in Tab. HI we perform a comparison to 
existing results in the literature. We compare our FRG 
results to those from proper-time flows (PTRG) (5t| and 
predictions from two different Monte Carlo simulations 
(MC) [13,|6l|. While [H provides for a statistical error, 
an estimate of the systematic errors of each calculation 
is not available. For instance, the PTRG results are on 
the one hand obtained from a larger truncation than the 
present one, but on the other hand do not rely on an ex- 
act RG flow [62| . The Nf = 4 calculation of 61 1 provides 
data for T cr /mf but no concrete error information. The 
Nf = 4 calculation ,60j provides data with a statistical 
error for T CT /ao which can correspond to T CI /mf within 
the conventions of [6^1 provided a multiplicative renor- 
malization is close to unity. Our results differ from MC 
[13 and PTRG results by 6% - 8%, which is a devia- 
tion slightly larger than our expected systematic error, 
whereas we find consistency with the MC data of 
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N f = 4 


JV f = 4 


AT f = 12 


AT f = 12 


AT f = 12 




(MC) m 


(FRG) 


(MC) 


(PTRG) m 


(FRG) 


T„(iV f ) 

7?lf 


0.49 


0.52 


0.66± 0.05 


0.57 


0.62 



Table I: Comparison of our results (FRG) for T CT (N{) to those of proper-time renormalization group flows (PTRG) [k| and 
Monte Carlo simulations (MC) [oOl. l6lT|. The conversion of the lattice data of [(5(| for T CI /a for ATf = 12 to our observable 
T cr /mt assumes that a potential multiplicative renormalization constant is close to unity. 




Figure 3: The phase boundary obtained from our flow equa- 
tions (red triangles corresponding to set I and orange tri- 
angles corresponding to set II, respectively) in comparison 
to available results from proper-time renormalization group 
(PTRG, Nf = 12) [H and Monte Carlo simulations (MC, 
Nf = 4, 12) [10, [6l[. 



within error bars. On the other hand, our curve appears 
to be a compromise between the existing results so far, 
sec Fig. El 

Let us now quantify the phenomon of local order intro- 
duced in the preceding discussion. As mentioned already, 
even above the critical temperature there is a tempera- 
ture regime T CI (N { ) < T < T*(Nf), where the flow first 
enters the broken regime at intermediate scales k but 
then re-enters the symmetric regime in the limit k — > 0. 
This implies that the flow exhibits a (nonunivcrsal) finite 
condensate at these intermediate scales which is even- 
tually depleted by the long-range bosonic fluctuations. 
Identifying 1/fc with a typical length scale, this interme- 
diate condensate can be interpreted as the existence of 
local or short-range order on these typical length scales in 
the system. 4 Only for higher temperatures T > T*(N{), 
no such local order is observed and the system is in the 
disordered state on all scales. Our results for T*(Nf) 
from numerically solving the flow are depicted in Fig. 0] 
(dashed line) in comparison with the critical temperature 



4 For example, the inverse IR regulator scale 1/fc can be roughly 
associated with the side length L of the spatial volume in lattice 
simulations [63ll , 



Figure 4: Local order in the Gross-Neveu model: the shaded 
local-order region is bounded by the critical temperature 
T CT (Nf) (solid line) and the temperature T*(Nf) (dashed line) 
above which the system is in the disordered regime on all 
scales. In the local-order region, a non-vanishing condensate 
is found at intermediate k scales, which is ultimately driven 
to zero by the long-range bosonic fluctuations. Interestingly, 
the ATf = 1 Gross-Neveu model shows an enlarged tempera- 
ture window compared to theories with ATf > 1. 



T CY (Nf) (solid line). It is interesting to observe that the 
Nf — 1 Gross-Neveu model shows an enlarged temper- 
ature window for local order compared to theories with 
Nf > 1, see Fig. [4] This might be due to the fact that 
already the UV fixed point potential is closer to the bro- 
ken regime. The temperature T*(Nf) lies always below 
the mean field critical temperature (dotted line). This 
is only natural, as the local-order phenomenon arises 
from an interplay of fcrmionic and bosonic fluctuations 
whereas the mean-field approximation ignores bosonic 
fluctuations completely, see our discussion above. 

At this point, we would like to add that a RG study 
of the Gross-Neveu model in a widely-used local poten- 
tial approximation (i. e. setting Z a = 1 and = 1) is 
not meaningful since it is spoilt by a dramatic parameter 
dependence. To be more specific, the parameters can be 
adjusted such that we obtain the same fermion mass in 
the vacuum for all Nf but the well-known exact result for 
the critical temperature in the large- Nf limit is neverthe- 
less not reproduced correctly, i. e. it differs, e. g., by a 
factor of two from the well-known result. This observa- 
tion can be traced back to the fact that the fixed-point 
structure of the theory is oversimplified in this approxi- 
mation. 
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We now turn to potential non-analyticities of the large- 
Nf limit. On the one hand, the large- Nf limit discussed 
in Sect. Qll]per constructionem predicts that the finite- 
temperature phase transition exhibits a mean-field be- 
havior with mean-field critical exponents. On the other 
hand, general universality arguments suggest that at any 
finite Nf, the phase transition should be in the 2d-Ising 
universality class, since the scalar order parameter mea- 
sures the breaking of a discrete Z2 symmetry. Also, 
close to the phase transition, the standard dimensional- 
reduction mechanism of the Matsubara formalism should 
be at work, distinguishing the bosonic zero-mode as the 
driving mode for the approach to criticality. 

The RG flow has to be able to resolve the obvious ten- 
sion between these two apparently contradictory view- 
points. One possibility is that this contradiction is re- 
solved by a non-analytic dependence of the critical phe- 
nomena on N{ potentially arising from a failure of the 
suppression of IR fluctuations with increasing Nf [6^ ]. 
The other, less dramatic possibility is that the size of the 
critical region close to the phase transition scales (pos- 
sibly non-analytically) with increasing Nf to zero, leav- 
ing us with the mean-field predictions in the Nf — > 00 
limit. Evidence for the latter possibility has been found 
in [35[ and [36[. In (35|, 2d-Ising scaling has been ver- 
ified by measuring the critical exponents very close to 
the phase transition with lattice simulations for Nf = 4, 
12, 24. On the other hand, beyond this regime a region 
obeying mean-field scaling has been identified and the lo- 
cation of the break-down of mean-field scaling in terms 
of a value for the condensate ctmf nas been determined. 
More precisely, for temperatures T < T cr correspond- 
ing to a condensate near <jq > cto.mFj mean-field scaling 
is observed, whereas strong deviations from mean-field 
scaling and a crossover towards the 2d-Ising regime oc- 
curs for a < ci m f ; c orresponding to temperatures closer 
to T cr . In Ref. [35|], the value of ctmf has been shown 
to scale as ctmf ~ 1/N^ MF . The numerical value for the 
exponent Xmf — 0.51 has been measured, nicely corre- 
sponding to an analytical estimate xmf = 1/2 which was 
derived from a simple estimate for the validity region of 
the mean- field approximation [35] ]. 

In [3(| , generalized large- Nf techniques have been used 
to compute the scaling of a fiducial relative temperature 
T — T cr (Nf) inside the Ginzburg critical region as a func- 
tion of Nf, yielding T - T CI (N { ) - 1/N? G with x G = 1. 
It was further argued that a scaling of the mean-field re- 
gion measured in terms of ctmf with xmf = 1/2 and of 
the Ginzburg region measured in terms of the tempera- 
ture with xq = 1 is compatible with each other within 
large- Nf considerations. 

In the present work, we go one step further and com- 
pute the condensate for a wider range of temperatures, 
resolving both the mean-field as well as the Ginzburg 
region. In particular, we directly determine the size of 
the Ginzburg region around the phase transition. We 
do so by quantifying the regime where 2d-Ising scaling 
is observed. For consistency reasons, the Ginzburg re- 
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Figure 5: The normalized order parameter squared as 
a function of temperature close to the phase transition. 
From top to bottom, the curves correspond to Nf = 
1, 2, 4, 6, 8, 10, 20, 30 , 100. For small Nf, the data is not well 
described by mean-field theory corresponding to a linear func- 
tion in this plot. However, as the flavor number increases, the 
temperature dependence is predominantly linear, as expected 
from the mean-field scaling with /3mf = 1/2. Deviations to 
linear scaling occur only upon approaching the critical tem- 
perature. As in [35|] , this suggests an increase of the width of 
the mean-field scaling region with increasing flavor number. 

gion must scale to zero with Nf 00 as fast as or even 
faster than the validity region of the mean-field approxi- 
mation. The size and scaling of the Ginzburg region is an 
important piece of information for experimental searches 
for the true critical region of the 2nd order phase tran- 
sition. In order to identify the Ginzburg region, we use 
the following strategy: The chiral condensate near the 
second-order phase transition scales as 

(a)~\T-T CI f, (31) 

where (3 in the exact solution should correspond to the 
exact 2c?-Ising critical exponent f3 = 1/8. The present 
FRG study indeed boils down to that of a 2c? scalar Z2- 
symmetric model near T cr due to dimensional reduction, 
guaranteeing that we are analyzing the correct univer- 
sality class. However, our simple truncation is not capa- 
ble of reproducing the exact Onsager critical exponents. 
Therefore, we identify the Ginzburg region by identi- 
fying the region in the phase diagram where Eq. (|3~lj) 
holds with a value for /3 which matches that of a 2c? 
Z2 theory in a cr 4 -truncation at NLO in a derivative 
expansion (with a linear regulator). For this, we find 
Arunc = 0.0931 using completely standard functional 
FRG techniques [55|. At fixed Nf, a one parameter fit 
of the form (a) = A\T - T cr |' 3tru » c with fit parameter A 
is then applied to a pre-selected region very close to the 
critical temperature, where the power-law is satisfied by 
the FRG data. For some flavor numbers, the FRG data 
is not in perfect agreement with the power-law, possi- 
bly due to inaccuracy in the determination of the critical 
temperature. Our predicted values for the size of the 
Ginzburg region should thus be taken with some care. 
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Figure 6: Double-log-plot of the width wp{N t ) (cf. Eq. (^21 1 
of the Ginzburg region close to the finite temperature phase 
transition as a function of the flavor number Nf. The tolerance 
criterion increases from bottom to top. The lowest data set 
corresponds to 2%, while the topmost corresponds to the 6% 
criterion. Within the Ginzburg region, 2d Ising scaling can 
be observed. The Ginzburg region is suppressed for larger 
flavor number with a powerlaw ~ 1/N^ G with varying from 
xq ~ 0.93 to ig — 1-06 (blue lines). 



Next, we define the temperature Tq(N{) characterizing 
the size of the Ginzburg region by requiring the relative 
deviation of the one-parameter fit to the FRG data to be 
less than a few percent inside that region. More precisely, 
we performed the determination of Tq(N{) from our data 
with several criteria corresponding to 2%, 3%, 4%, 5% 
and finally 6% tolerance. 

Repeating the analysis for various Nf, we find that we 
can fit the corresponding width wp (Nf) of the Ginzburg 
region scales to a powerlaw 



wp(Nf) 



T G (N{) - T CI {N{) 



1 



(32) 



where the exponent obtained from the fit varies from 
xq ~ 0.93 to xq — 1-06 in the range of tolerances from 
2% to 6%. This apparent definition dependence of xq 
may also originate from the influence of subleading pow- 
ers of 1/Nf which need not be vanishingly small in the 
fit regime. Still, our result is well compatible with the 
generalized large- Nf result xq = 1 [36j - In this sense, it 
is remarkable that this scaling law appears to hold down 
to rather small values of Nf. A plot of the scaling laws 
is shown in Fig. [6] Technically, the width becomes more 
difficult to determine for increasing Nf, as the resolu- 
tion along the temperature axis has to be finely adapted. 
In Fig. [6] we show our data in a range from Nf = 1 to 
Nf = 30. 

Comparing our result for the scaling of the Ginzburg 
region, xq ~ 1 with that of [35[ for the scaling of the 
condensate value at the boundary of the mean-field re- 
gion xmf — 0.51, the two different values can be con- 
nected by a plausibility argument: our criterion derives 



Figure 7: Plot of the classical regime (dashed black) below the 
phase boundary in the many-flavor phase diagram. The criti- 
cal region begins where the renormalized thermal mass of the 
bosonic fluctuations equals temperature, T cl = m(T cl ). The 
physics in this region is dominated by the bosonic Matsubara 
zero mode. 



from a quantity scaling with temperature, cf 
whereas 



Eq. 432), 

studied the scaling of the condensate in the 
mean-field region. Since mean-field theory predicts a re- 
lation between both with an exponent Pmf = 1/2, it is 
reassuring that the difference between J35| for the scal- 
ing exponent of the condensate value at the boundary of 
the mean-field region xmf and the scaling exponent of 
the Ginzburg region xg are separated by an amount on 
the order of /3mf- A formal large- Nf argument (3(| also 
points into the same direction. 

Let us conclude this part with a word of caution: Even 
though the properties within the critical Ginzburg re- 
gion are universal, the boundary of this region measured 
in terms of Tq as well as the boundary of the mean- 
field region may not exhibit the same degree of univer- 
sality. This is, because Tq is determined by the tempera- 
ture where the shorter-range thermal fermion and higher- 
Matsubara boson fluctuations kick in. Such sources of 
non-universality should even be more pronounced for the 
location of the boundary of the mean-field region and the 
corresponding exponent xmf, since the mean- field region 
is even further apart from the universal critical regime. In 
fact, we observe an even stronger definition dependence 
of our results for imf which nevertheless is roughly of 
the same size as the lattice value (35j | . 

Finally, let us analyse the size of the so-called classical 
regime, where the dimensional reduction, i.e. the dom- 
inance of the bosonic zero mode sets in in the broken 
phase. Following (35|, we define the onset of the clas- 
sical regime by the value T cl of the temperature where 
this temperature equals the renormalized bosonic mass 
T cl = m(T cl ). In the region T cl < T < T cr , the bosonic 
Matsubara zero mode contributes significantly to the 
flow. This classical region is shown as a shaded area 
in Fig. [7J Its width appears to be approximately inde- 
pendent of the flavor number. More quantitatively, we 



11 



find that the width w c \(Nf) of this region, measured in 
units of the zero-temperature fermion mass, is 

r cr (jv f )-r cl (jv f ) 

Wci(A'f) = — const. , (33) 

mf 

aquires values in the range w c i — 0.15... 0.25 with the 
smaller values holding for Nf < 10. Therefore, the classi- 
cal regime remains finite even for Nf — > oo. This simply 
reflects the fact that the massless bosonic Matsubara zero 
mode will always exist near the phase transition even in 
the presence of an arbitrarily large but finite number of 
(thermally massive) fermion modes. This phenomenon 
remains completely invisible in the strict large- Nf limit. 
This observation indicates that the formal large- Nf limit 
has to be considered with some care: taking the limit 
Nf — s- oo on the level of the path integral does not appear 
to be equivalent for all physical observables to solving the 
theory including bosonic fluctuations and then taking the 
limit Nf — > oo. 

VI. CONCLUSIONS 

We have investigated the many-flavor phase diagram 
of the (2 + l)-dimcnsional Gross- Neveu model at finite 
temperature. At zero-temperature, the model undergoes 
a quantum phase transition as a function of the coupling 
strength, separating a massless from a massive phase. 
The critical coupling marks a non-Gaufiian fixed point 
which renders the model asymptotically safe, i.e. UV 
complete despite its perturbative non-renormalizability. 
As this fixed point has only one relevant direction, all 
long range observables can in principle be predicted once 
a single physical scale is fixed. Hence, also the many- 
flavor finite-temperature phase diagram for the massive 
model is also a pure prediction of the theory, if expressed, 
for instance, in units of the vacuum fermion mass mf. As 
mass generation as well as symmetry restoration at fi- 
nite temperature is a non-perturbative process, we have 
used the functional RG for a determination of various 
aspects of the phase diagram. In a conceptually com- 
paratively simple framework, the functional RG gives 
direct access to universal quantities in the critical re- 
gion near the phase transition as well as non-universal 
properties. In particular, our predictions for the finite- 
temperature many-flavor phase boundary can be tested 
by other non-perturbative approaches, such as lattice 
Monte-Carlo simulations. 

In addition to the phase boundary, we have particu- 
larly concentrated on the flavor-number dependence of 
the phase-diagram regions of (i) local order, (ii) the 
Ginzburg critical region and (iii) the classical region, all 
of which correspond to beyond-mean-field or finite-iVf 
phenomena. 

The region of local (or short-range) order at tempera- 
tures right above the phase transition characterizes tem- 
perature where a finite condensate can be found on inter- 
mediate length scales. In addition to possible influences 



on the spectral function, an understanding of this region 
is particularly important for finite- volume studies, since 
local order can easily be confused with long-range order 
for too small volumes. It would be indeed interesting to 
further analyze this question, e. g., by directly employing 
lattice Monte-Carlo simulations. 

The Ginzburg critical region corresponding to the 
range of temperatures where the phase transition shows 
universal 2d Ising scaling is especially relevant for an un- 
derstanding of the approach to the large- Nf limit. The 
non-analytic jump of the critical exponents from 2d Ising 
to mean field at Nf — > oo is realized by the size of the 
Ginzburg critical region vanishing with 1 /Nf° . Our re- 
sult xq ~ 1 is in accordance with other values in the 
literature (35l. l36j). 

By contrast, we find that the classical regime, where 
the dynamics of the bosonic Matsubara zero-mode sets 
in below T cr , remains finite at any Nf. The existence of 
this mode cannot be inhibited by a numerical dominance 
of a large number of flavors. The formal large- Nf limit 
should therefore be viewed as a pure projection technique 
onto the fermionic fluctuations rather than a means for 
quantifying the dynamical prevalence of fluctuations. 

With respect to the applications of the Gross-Neveu 
model in condensed-matter systems, the Nf = 2 model 
is phenomcnologically most important. Whereas some 
quantities such as the many-flavor phase boundary may 
semi-quantitatively be estimated by higher-order stud- 
ies in a large-TVf expansion even down to smaller val- 
ues of Nf, others such as the region of local order show 
strong variations particularly near Nf = 2. This inhibits 
a general statement about the in-principlc applicability 
or inapplicability of large- Nf expansions for physical phe- 
nomena such as instabilities in graphite or graphene, or 
secondary d- to d + is-wave pairing transitions in nodal 
cf-wave superconductors. 
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Appendix A: Threshold functions 

The regulator dependence of the flow equations is en- 
coded in dimensionlcss threshold functions which are as- 
sociated with 1PI diagrams incorporating bosonic and 
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fcrmionic fields. In this work we have employed a covari- 
ant linear regulator [57j for which the threshold functions 
at zero temperature are known in closed form, see, e.g., 
Ref. [H. 

At finite temperature, the loop momentum integrals 
have to be performed by replacing the time-like mo- 
mentum component by the corresponding fcrmionic or 
bosonic Matsubara frequencies. The resulting threshold 
functions thus acquire a dependence on the dimensionlcss 



temperature parameter r = 2irT '/k. For the treatment of 
loop diagrams with mixed fermionic and bosonic propa- 
gators in the loop, we follow the prescription of |66[ in 
order to achieve consistency with the derivative expan- 
sion. With the same techniques and definitions as for 
T = and with v^ 1 = 2 d+1 n d ^T(d/2) accounting for 
Euclidean volume factors, we obtain the following finite 
temperature threshold functions 



l n {T,UB;V<r) = o 7T / \\ TC B,n) + V) {Ga {TC B ,n , V, U B ) ) =; ; 

2 Vd 27r nek-* Za > k 

V T f°° d 3 

Vd 2n ^ J ° 

(l + r F )d t {Z atk r F ) 



Z 



li F m( T ,u F ,ujB;iU,ifc) = ~7:~~ — ~7T~ / d yy^d t {(Gi,(Tc F .n,y-,uj F )) n (g (7 (Tc F . n ,y;uj B )y n } 7 

Z Vd Z7T ' — ' In 

T f°° d 1 ( \ ^ 



{F)d, \ 



n£2 

1 Vd-i d t 

2 v d d- 



(F)d, s lwd-i d T 



d t poo d i f d 

JZ7l2^Yl / d y y^ i( TC F,n) 2 + y)d t ( Q-[(l + r F )g^(Tc F , n ,y;u} F )\ 



2 v d 



i/2 



-Q^Ga(TCB,n,y,U B ,2) 



m[ F 2 B)d {T,u F ,u B ; Ti^,T] a ) = ~\-^j—J^. / dyy^'dt |(1 + r F ){G^,(y;u F )) ( -^Gaiy^B) 



dy~ 



where the operator dt is defined as 



St\ a = [-2((rc„) 2 + y y B -. qa r B \^-, (Al) 

dtU = [-2(( T c n ) 2 + y)r' F - ^rp] (A2) 

depending on whether it acts on bosonic or fcrmionic propagators. Here, the Matsubara frequencies are encoded by 
the variable c n , acquiring the values c n = n for bosonic terms and c„ = n + ^ for fermionic terms. "Covariant" 
optimized shape functions now read 

r B ((rc n ) 2 +y) = ( ^ - - - ij 6(1 - [(rc») 2 + y]), (A3) 



r F {{rc n f+y) = ^ ^—^ - - } 1 j 9(1 - [{rc n f + y]). (A4) 

The threshold functions for the linear regulator can be factorized in terms of thermal and vacuum contributions. In 
the following, full thermal threshold functions and their vacuum counterparts arc distinguished by the dependence on 
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the temperature parameter r. We find 

IHt^b;^) = (^(t)- ??ct ^(t))/> b ;0), 



li n F nA T : W F , W B ; 7fy , 7^ ) 



1 + CJ F / \1 + Wfi 

fnf )d (T,uj F ;ri^) = 4 F ' d (T)m\ F ' d (w F ;0), 



1 + ojf 



1 - w 



*4(r) 



1 



1 + 

2 



1 ~~ . / \ . 1 (F)d, s 

-J-kW + ^o (r) ;vi + wi? 

m lm( T i w B,l ! WB,2;i,r) = Sq ( T ) m n,m ( W B : 1 , W B , 2 ; 0) , 

m™^,^;^,^) = (>(t) - ^^(r)) mg B)d (cj F , o; B ; 0, 0). 
The thermal threshold factors are themselves given by the following Matsubara sums 

4(r 
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1 2tt 



d— 1 
-2\ S~ 



d-l 

~2~ 



d-l 



d+1 
r 2\ 2 



1 - |cf,„t| 2 Fi -1/2, 



d-l d+l c| „t 2 - 1 

~ ; ~ ; 



if £ W - (1 - cV 2 ) ^ Ic^|- 2 Fr f 1/2, 1±±; f*±* ; , 



with the hypergeometric function 2-F\(a, 6; c; z) [6f 



In the following, we give closed expressions for the ther- 
mal threshold factors, specifically for three spacetime di- 
mensions (d = 3), which can be obtained by applying the 
Poisson rcsummation formula 



£/(») = £ / <*<?/(<?) e" 2 ^, 



(A5) 



and the identity [67[ 

Li„(-z)-(-l)«- 1 Li„(-l/z) = - ( ^B n (^4 + |) , 

(A6) 

for 2 ^ (0, 1) where B n (z) are Bernoulli polynomials and 
Li n (z) denotes the nth polylogarithm. With 



In 



(V 27ri/r ) = -27ri + 27ris F (r), 



In [- C - 2 ^ T ^ 



27ri „ , , , , . _. 
i7T + 2itisb(t), (A7) 

T 
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Figure 8: Comparing j-Si (^) (blue) and the approximant 

X]fc=o f {Jk+1/2 (^t)) 2 (magenta) as a function of r for the 
Poisson indices I = 1,2,3,4. The Poisson resummation over 
I can be done in closed form for the approximant. The de- 
viation for t — > is immaterial since the sine integral enters 
the thermal threshold factor with a prefactor r and the flow 
is stopped at finite r. 



4(r) 
s d (r) 



(r) 

4 F) V) 

s d (r) 
s d (r) 

U(t) 



1 i (-(r + 2ts b (t)) (-3 + r 2 s B (r) + r 2 s B (r) 2 )) 



d — 1 I'd 2n 3 
d 1 1 



-(1 + 2* B (r)) (l5 - t 2 (10 + r 2 ) s b (t) + 



d 2 - 1 w d 2tt 15 

2r 2 (-5 + r 2 ) SB (r) 2 + 6t 4 s s (t) 3 + 3r 4 Ss (r) 4 ) , 



1 ir.s F (r)(l2+(l + 12)r 2 -4r 2 Si ,(r) 2 ), 



d — 1 Wrf 27r 6 
d 1 1 



4s F (r) + (4 + r 2 ) s f (t) - t 3 s f {t) 2 - At 2 s f {t) 



d-l Vd 2?r 12 
2t 3 s f (t) 4 + r 2 (s F (r) - ts f {t) 2 - 4s F (r) 3 + 2ts f (t) 4 )^, 

d 1 1 ( TSF ( T )(24 + 40r 2 + 7r 4 - 40r 2 (4 + t 2 )s f t 2 + 48t 4 s f t 4 )) 

cr — 1 27T 12U 

1 ^ts f (t) (4 + r 2 + t 2 s f (t)(-t + 2s f (t)(-2 + ts f (t)))) , 



d 2 - 1 u d 2?r 3 
d Ud-l 1/1 



ts f (t)(-12 - r 2 + 4t 2 s f (t) 2 ) + -t(32s f (t) + 2r(-12 - r 2 )s F (r) 2 
ci — 1 Vd Ztt y b 3 

1 



4r 3 s F (r) 4 - — r(-32s F (r) + 2r(40 + 20r 2 + 7t 4 )s f (t) 2 
40r 3 (-2 - r 2 )s F (r) 4 + 32t 5 s f (t) 6 )) + . . . ^ . 



The expression for <4(t) is only approximate since the cally. However, after performing the Poisson integra- 
Poisson resummation could not be performed analyti- tion, the resulting epxression contains a sine integral, 
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which can be expanded in terms of Bcsscl functions, 

Si(z) = T^Y^k=o {Jk+i/2 (§)) ■ Above, we expanded to 
k = 2 and cross-checked with numerical evaulation of 
t±{r) : see also Fig. El The agreement is satisfactory. In 



fact, we have checked that the difference between our ap- 
proximate expression and the exact result does not lead 
to significant differences in our results for the phase dia- 
gram. 
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